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In this paper we study equations of magnetic hydrodynamics with a stress tensor. We interpret 
this system as the generalized Euler equation associated with an abelian extension of the Lie 
algebra of vector fields with a non-trivial 2-cocycle. We use the Lie algebra approach to prove 
the energy conservation law and the conservation of cross-helicity. 



1. Introduction. 

Magnetic hydrodynamics (MHD) describes evolution of a fluid or plasma, carrying 
a magnetic field. This theory is used to model the processes in the Solar corona [17], as 
well as to design tokamaks [10]. There are numerous books treating various aspects of the 
subject, see e.g., [8], [15]. 

The MHD equations are derived from the Euler equation of motion of an incompress- 
ible fluid and the Maxwell's electrodynamics equations, and describe evolution of a fluid 
with the velocity vector field v and the magnetic field B: 

^ = _( v • V)v + (B • V)B - Vp, 

^ = -{v,B}, (1.1) 
div(v) = 0, div(B) = 0. 

In this paper we study another system of PDEs, where we add extra terms into the 
evolution equation for the velocity: 

^ = - (v .V)v + (B.V)B + E ||vg-V P . (1.2) 

The additional terms here are of the third degree in derivatives and one can draw certain 
parallels between (1.2) and the Korteweg - de Vries equation 

u t = uu x + u xxx . (1.3) 

We interpret the extra terms in (1.2) as a contribution of a stress tensor 

k% dxj dxk ^ ^ 

j J 

This stress tensor is not symmetric, which indicates that the particles of the fluid should 
possess electric or magnetic momentum. 
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In spite of the additional higher-order terms, the new system of PDEs retains many 
features of the original system. In particular, we show that it still admits the Alfven wave 
solutions. 

We follow an approach developed by Arnold to give an interpretation of the MHD 
equations with the stress tensor as a generalized Euler equation. A generalized Euler 
equation is an equation for the geodesies on a (possibly infinite-dimensional) Lie group 
supplied with a Riemannian structure. This equation describes the evolution of a tangent 
vector of the geodesic in the Lie algebra of the Lie group (see Section 3 for details). 

The Lie algebra that corresponds to the MHD equations with the stress tensor is an 
abelian extension g(r) of the Lie algebra of the divergence zero vector fields twisted with a 
non-trivial 2-cocycle r. This Lie algebra was studied in the framework of the representation 
theory of the toroidal Lie algebras and the cocycle plays a prominent role there. 

Infinite-dimensional groups associated with abelian extensions of the Lie algebra of 
vector fields are discussed in [7]. 

The Lie algebra q(t) has nice properties and this translates into nice properties of 
the PDEs. In particular we establish the energy conservation law and the cross-helicity 
conservation for MHD with the stress tensor. 

Ovsienko and Khesin showed in [16] that the generalized Euler equation for the Lie 
algebra of vector fields on a circle yields the non-linear wave equation ut = uu x , while 
incorporation of the Virasoro cocycle into the Lie algebra leads to the Korteweg - de 
Vries equation (1.3). In a way, the present paper may be viewed as a higher-dimensional 
generalization of [16]. 

The paper is organized as follows: in Section 2 we discuss the properties of the system 
(1.2), derive the expression for the stress tensor, list the conservation laws and describe 
the Alfv§n wave solutions. In Section 3 we review the generalized Euler equation for an 
arbitrary Lie algebra and we apply this method in Section 4 to an abelian extension of 
the Lie algebra of the divergence zero vector fields, deriving (1.2), and establishing the 
conservation laws in a purely algebraic way. 

2. Magnetic hydrodynamics with a stress tensor and its properties. 

Evolution of an incompressible fluid carrying a magnetic filed is given by the equations 
of magnetic hydrodynamics (MHD): 



Here B is the magnetic field, v is the velocity vector field of the fluid and p (pressure) is an 
auxiliary function which is chosen in such a way that the equation div(v) = is satisfied. 
The formal dot product v • V represents the differential operator 




div(v) = 0, div(B) = 0. 



(v • V)v + (B • V)B - Vp, 



(2.1) 
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The Poisson (Lie) bracket {v, B} of two vector fields v = ^ Vj (x) and B = ^ Bj(x) 

3 3 J 



is given by 

{v, B} = ad(v)B = (v ■ V)B - (B • V)v = £ (v^ - . 

In the three-dimensional space, the first equation in (2.1) may be also written as 

<9v 

— = v x curl(v) + curl(B) x B — Vp. 

at 

For a conducting medium, the curl of the magnetic field is equal to the electric current 
(Ampere's law), and so the expression curl(B) x B represents the Lorentz force, with which 
the magnetic field acts on the current. 

The configuration space for a flow of an incompressible fluid is the group of volume- 
preserving diffeomorphisms. In his remarkable paper [2], Arnold interpreted the Euler 
equation for an ideal fluid as the geodesic equation on this Lie group. The geodesic 
equation describes the evolution of the tangent vector of the geodesic curve, and this 
tangent vector belongs to the Lie algebra, which is the Lie algebra SVect of the divergence 
zero vector fields, in case of the group of volume preserving diffeomorphisms. 

Using Arnold's method, Vishik and Dolzhanskii [19] (see also [12]) showed that the 
MHD equations (2.1) also may be interpreted as a geodesic equation for a certain infinite- 
dimensional Lie group. The Lie algebra that is used to write this equation is the semidirect 
product 

SVect ©fi 7^° 

of the Lie algebra of divergence zero vector fields SVect with its dual space - the factor 
O 1 /cZfi of the differential 1-forms modulo the exact 1-forms. We review this construction 
in detail in Section 3. 

The Lie algebra that is associated with the MHD equations, 

g = SVect © tt^dtt , 

has recently attracted much interest in representation theory. It turns out that represen- 
tations of g may be used for constructing modules for toroidal Lie algebras (see e.g., [4], 
[5], [9]). It was also discovered that this Lie algebra has an important deformation - the 
Lie bracket in q may be twisted with a Virasoro-like 2-cocycle r. The twisted Lie algebra 
q(t) still has nice properties, and its representation theory is even better than of g itself. 

In this paper we study the system of PDEs that comes from the geodesic equation for 
fl(r). In Section 4 we show that this Lie algebra yields the following system of PDEs: 

^ = -(v.V)v + (B.V)B + E||V§|-Vp, 

f = -{v,B}, l " ( 2 - 2 ) 

div(v) = 0, div(B) = 0. 
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Let us discuss the properties of this system of PDEs. 
First of all, we note that the new term 



^ = £^(^1 (2-3) 

can also be written as 

Edvj^dB; 



since the difference of the two expressions is the gradient of ^ and may be absorbed 

i,j 1 3 

into Vp. 

It is curious to note here that in the one-dimensional case the passage from the Lie 
algebra of vector fields to the Virasoro algebra leads to the transition from the non-linear 
wave equation 

u t = uu x 

to the Korteweg-de Vries equation 

Ut = uu x -\- U XXXl 

as shown by Ovsienko and Khesin in [16]. 

Just like the dispersion term u xxx in the KdV, the new term ff 1 ^"^" that we get 

i,j 3 

in (2.2) has a triple derivative in x. 

Next we are going to show that (2.2) describes magnetic hydrodynamics with a stress 
tensor. Indeed, for a stress tensor Tm, the equations on the velocity field in (2.1) will 
become (see e.g. Section 1.7 in [18]) 

— = -(v • V)v k + (B ■ V)B k + £ (2.4) 

i 

Proposition 1. The system (2.2) describes magnetic hydrodynamics with a stress 
tensor. The stress tensor T ki may be written as 

Q x Q Xk ft.5i) 

3 

or as 

3 J 

or as a linear combination aT^ + (3T' ki with a + (3 = 1. 
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Proof. We will prove the statement of the Proposition for tensor T ki given by (2.5). We 
write the contribution in (2.4) from the stress tensor (2.5): 



Since div(B) = 0, the first term in the right hand side vanishes, and we get precisely the 
first equation from (2.2). 

The stress tensor (2.5') will yield the additional term written in the form (2.3'). The 
proof in this case is completely analogous. 

We point out that the stress tensors we obtain here are not symmetric: T ki ^ T ik . 
Asymmetric stress tensors occur when the particles of the fluid are polar, i.e., possess 
electric or magnetic momentum (see Chapter 8 in [18]). The derivation of the stress tensor 
from the first principles is rather delicate (see e.g., Chapters 7 in [18]), and I am unable 
to describe precisely the physical situations when the stress tensors (2.5) or (2.5') would 
occur. 

Introduction of the term (2.3) which has the 3rd order in derivatives, into the equations 
will clearly change the behaviour of the solutions in a substantial way. I view of that, it 
is quite surprising that the conservation laws of magnetic hydrodynamics still hold for the 
system (2.2). 

Our next goal is to study the conservation laws for magnetic hydrodynamics with the 
stress tensor (2.2). However before we do that, let us discuss the class of solutions that we 
consider here. 

We require that the functions Vi(x), Bi(x) are defined in a domain D C M n and belong 
to the intersection of the Sobolev spaces n Hq(D). We recall that the space Hq(D) is the 



closure in the Sobolev space H k (D) of the functions of class C k with compact support (see 
e.g. [1]). This will ensure that all functions and their partial derivatives of all orders are 
square-integrable (belong to L 2 (D)) and satisfy the vanishing conditions on the boundary 
of D (if D is unbounded this means that at infinity the functions go to zero faster than 
the inverse of any polynomial). Choosing this class of functions will allow us to carry out 
integration by parts with the boundary term vanishing. 

Alternatively, we may consider periodic boundary conditions, as it is often done in 
turbulence theory. 

There are several conserved quantities for the MHD system - mass, momentum, mag- 
netic helicity, as well as energy and cross-helicity. It turns out that all of these are also 
conserved for the MHD with the stress tensor. Since we consider the case of an incom- 
pressible fluid, the conservation of mass (volume) holds trivially. The derivation of the 
conservation of magnetic helicity involves only the evolution equation on magnetic field B, 
and is exactly the same for both systems. Let us prove the conservation of momentum for 
the new system. 

Proposition 2. The total momentum is a conserved quantity for the MHD system 
with the stress tensor (2.2): 



E 



dT ki 



dxi 
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Proof. With the help of Proposition 1, we can write the first equation using the 
divergence operator: 



dt dx % 

i 

By the divergence theorem, we get 



= ^2 ~7T~. (~ v i v k + B i B k + T ki ~ pSki) ■ 



^- [ v k {x)dV = (f R ki ei ■ dn, 
ot J D Jq D 



where R is a 2-tensor 

Rki = —ViVk + BiB k + T ki — pSki, 

n denotes the unit outward normal vector, and e^'s are the standard basis vectors. Since 
the vector fields we consider vanish on the boundary of D, the last integral is zero. 

Next we state the corresponding theorem for the conservation of energy and the cross- 
helicity conservation: 

Theorem 3. The system (2.2) of magnetic hydrodynamics with stress tensor (2.5) 
in a domain DcK" with appropriate boundary conditions (see discussion above) has the 
following two first integrals: 

/ Vj(x) 2 + Bj(x) 2 dV = Const (energy conservation) (2.6) 



and 



/ Vi(x)Bi(x)dV = Const (cross-helicity conservation). (2.7) 



We will give the proof of this theorem in Section 4. These conservation laws will be 
derived from the properties of the Lie algebra $j(t). 

It is interesting to note that unlike the case of the Navier-Stokes equation, introduction 
of the stress tensor in (2.2) does not lead to the dissipation of energy, and the energy 
conservation law still holds. 

For the topological interpretation of helicity and cross-helicity, see [13], [14]. 

In conclusion of this section we are going to show that system (2.2) admits Alfv§n 
wave solutions. 

Alfven waves solutions are obtained as a perturbation of a steady-state constant so- 
lution v(x) = 0, B(x) = B . If we take an expansion v = v(x), B = B + B(x) near this 
equilibrium state, we will get the following system: 



§ = -(v.V)v + (B -V)B + (B.V)B + £ftvg-Vp, 

§ = -(v • V)B + (B • V)v + (B • V) v, (2 - 8) 
div(v) = 0, div(B) = 0. 
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Next we set B = v. In this case the term Y ff^^l^r ls a g ra *dient of Y an d we 

i,3 3 % i,3 3 

can eliminate it by setting 

i,3 J 

Now the first two equations in (2.8) reduce to a single equation 

^ = (B -V)v. (2.10) 

Finally, by taking an arbitrary divergence zero vector field w(x), we can construct a 
solution of (2.8): 

v(x, t) = w(x + B i). 

This traveling wave solution is called the Alfv§n wave. The only difference with the classical 
MHD system (2.1) is the change of the pressure term (2.9). 



3. Generalized Euler equation. 

In this section we are going to review the geodesic equation approach to hydrodynamics 
developed by Arnold. In the key paper [2], Arnold gave an interpretation of the Euler 
equation for an incompressible ideal fluid 

/fi = -(v-V)v-Vp , , 

\div(v) = 0, {6A} 

from the perspective of infinite-dimensional Lie groups. He showed that the Euler equation 
may be interpreted as the geodesic equation on the group of diffeomorphisms, where the 
Riemannian structure on the group is given by the energy functional. 
We will describe this approach here following the book [3]. 

Let G be a Lie group (possibly infinite-dimensional), and let g be its Lie algebra. 
Consider a map from g to its dual 

^4 : 0^0* (inertia operator), 

such that it defines a positive-definite symmetric bilinear form on g: 

<X\Y >=±A(X)Y+±A(Y)X, X,Y e g. 

The corresponding quadratic form < X\X >= A(X)X is called the energy functional. 
The Lie algebra g acts on its dual space via the coadjoint action: 

(ad*(X)u) (Y) = -u([X,Y]) for X,Y G g, u G g\ 
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We assume that the space A(g) is invariant under the coadjoint action, and we make a 
convention that in what follows g* stands for A(g) (this is a slight abuse of notations 
since in the infinite-dimensional case A(g) is typically smaller than the formal dual of g). 
Since the kernel of A is trivial (otherwise, the quadratic form < X\X >= A(X)X is not 
positive-definite), then with the above convention the operator A is invertible. 

The bilinear form < -| - > may be left-translated from g = T e (G) to the tangent spaces 
at all points of G. This gives a Riemannian structure on G, and allows us to consider the 
geodesies on this group. 

Next we are going to write the equation for the geodesies on G, which describes the 
evolution of the tangent vector X e g to the geodesic curve. It turns out however, that it 
is easier to write the evolution equation for the covector u = A(X) e g* rather than for 
X itself. The generalized Euler equation is the evolution equation for u = A(X) which is 
written using the coadjoint action (see (6.4) in [3]): 

u t = -ad* (A -1 w) u. (3.2) 

When G = 50(3) this equation turns into the equations of motion of a rigid body 
with a fixed point. 

Let us discuss equation (3.2) in the context of fluid dynamics. Evolution of an 
incompressible fluid in domain D C IR n from time to time t is given by a volume- 
preserving diffeomorphism of D. Thus the group of the volume preserving diffeomorphisms 
G = SDiff(L>) is the configuration space for this motion. The Lie algebra of the group 
SDiff(D) is the Lie algebra of the divergence zero vector fields SVect(-D). As the energy 
functional we take the kinetic energy: 

(E^w^lE«.w^) = / D i:«.w a ^ (»■») 

The dual space for the divergence zero vector fields is the factor space of differential 
1-forms on D by exact 1-forms: 

SVect(D)* = n 1 (D)/dn°(D). 

The pairing between SVect(-D) and fi 1 (D) / <ifi° (D) is given by the integral: 

It is easy to check that exact 1-forms vanish on the divergence zero vector fields, and so 
the value of the integral on the right is independent of the choice of a representative in a 
class of 1-forms modulo dflP(D). 
The inertia operator 

A : SVect(D) -> Q 1 (D)/dfl°(D) 
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that corresponds to the energy functional (3.3) is written as follows: 



A ^2vi(x)-^j = ^v i {x)dx i 



It is possible to check that in these settings, equation (3.2) turns into the Euler equations 
of motion of an ideal fluid. 

The equations of the magnetic hydrodynamics (2.1) can also be obtained as a special 
case of the generalized Euler equation (3.2). To construct the corresponding Lie algebra 
we take the Lie algebra SVect(-D) together with its dual space: 

g = SVect(D) © Q 1 (D)/dtt°(D). 

The Lie bracket of two 1-forms is set to be zero. The Lie bracket of two vector fields is 
defined in the usual way: 



dXn 



■I -J 



dw-i d 



dvi d 



— Wn 



1 dxi dxj 3 dxj dxi ' 



(3.5) 



and the Lie bracket of a vector field with a 1-form is given by the Lie derivative: 



d 



dw 



(3.6) 



It is easy to see that the space dQ°(D) is invariant under the Lie derivative action, so the 
above formula may be taken modulo dfl°(D). 

The Lie algebra g is a semidirect product of SVect(-D) with its module O 1 (D) / <ifi° (D) . 
The space Q 1 (D)/dQ°(D) forms an abelian ideal in q. 

A really important feature of the Lie algebra q is the existence of an invariant sym- 
metric non-degenerate bilinear form. In contrast, the Lie algebra SVect does not possess 
such a form. The invariant form on g is defined as follows (cf. (3.4)): 



5><(*)bS- 



J2 w j( x ) dx 3 = I D Y, v i{ x ) w i{ x ) dV i 



Y / w j (x)dx j ] = 0. 



(3.7) 



One can verify that the bilinear form (3.7) satisfies the invar iance property: 

([X, Y]\Z) = (X\ [Y, Z]) for all X, Y, Z e 0. 

We can use this form to identify each element leg with a linear functional (X\-) 
in g*. It is well-known that when the bilinear form that is used to identify g* with g, 
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is invariant and non-degenerate, the coadjoint action becomes isomorphic to the adjoint 
action. In this case the generalized Euler equation takes form: 

X t = -ad(A~ 1 X)X, (3.8) 

where now X e Q and the inertia operator A now maps q to g. 

The generalized Euler equation (3.8) yields the equations of magnetic hydrodynamics 
(2.1) if we choose A to be the following involution on q: 



A (j2 v i( x )-£^j = H v i{ x ) dx h 



(3.9) 



Note that in the second equality we choose a (unique) representative in a class modulo 
dfl° satisfying = 0, so that the right hand side is a divergence zero vector field. 

i % 

We see that the inertia operator (3.9) satisfies A -1 = A, and the energy functional 
(AX\X) for X = J2 v i( x )~£r + S Wi(x)dxi E g is given by the integral 

i i 

[ $>(a0 2 + yV(a0 2 dV. (3.10) 

The Lie algebra 

= SVecteO 1 /^ 

appears in the study of toroidal Lie algebras ([4], [5], [6], [9], [11])- The representations 
of are an essential ingredient for constructing the representation theory of toroidal Lie 
algebras. It was discovered however that g admits a non-trivial deformation with a Q 1 /de- 
valued 2-cocycle on SVect, and one gets a better representation theory for the deformed 
algebra than for q itself. 

In the next Section we will describe this deformation of q and study the associated 
generalized Euler equation. 



4. Abelian extensions of the Lie algebra of vector fields. 

In the previous Section we have constructed a semidirect product q of the Lie algebra 
of divergence zero vector fields with O 1 /^ . It turns out that on the same vector space 

g = SVect ©O 1 /^ 

we may deform the Lie bracket in a non-trivial way. When we define the Lie bracket of two 
vector fields, we are going to add to the right hand side of (3.5) a correction term which 
has value in fl 1 (D)/dn°(D): 
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= E 



1,3 



dwj d dvi d 

V; —r^- — W 



dxj, dxj 3 dxj dxi 



In order to get a Lie bracket, r has to be a 2-cocycle on SVect(D) with values in 
Q 1 (D)/dQ°(D). 

The following cocycle plays an important role in the representation theory: 

\ 1 3 J / 1,3 J v ' 

This cocycle may be viewed as a higher-dimensional generalization of the Virasoro 
cocycle. Just as the Virasoro cocycle, it has a triple derivative in x, and in fact (4.2) 
reduces to the Virasoro cocycle for the Lie algebra of vector fields on a circle. 

We will denote the Lie algebra with the Lie bracket deformed by the cocycle r by 
fl(r). Note that 

( r ) = SVect(D) © Q, 1 (D) / dfl° (D) 

is no longer a semidirect product, but the subspace Q 1 (D)/dQ°(D) still forms an abelian 
ideal. As before, the action of SVect(D) on Q 1 (D) / dQ° (D) is given by the Lie derivative 
formula (3.6). 

Proposition 4. The bilinear form on jj(t) given by (3.7) is invariant. 
Proof. We need to establish the invariance property: 

{[X,Y]\Z) = [X\%Z]). 

There are three non-trivial cases to be considered: 

(i) X, Y e SVect(D), Z E O 1 (D) / dfl° (D) , 

(ii) X,Z E SVect(D), E G n 1 (D)/dfl°(D), 
(hi) X,Y,Z E SVect(D). 

We will verify the invariance only for the last case, since only this case will involve 
the cocycle r, and leave the first two cases as an exercise to the reader. 

Suppose * = £^^,E = y> J gf-,Z = y> fe g|-. Since (S Vect | S Vect) = 0, then 

i 3 k 

we get that ([X,Y] \Z) = (t(X,Y) \Z) and (X\ [Y, Z\) = (X\t(Y,Z)). 



We have 



([X,Y]\Z) = \^ du < d% "> 



^ dxj dxidx 




= I dui d2yj dV. 

Id , dxj dxidxk 

i,3,k 



Integrating by parts and using the fact that V] = 0, we get 



([X,Y]\Z)=-[ 



dw k duj dvj . 
dxi dxj dxk 
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On the other hand, 



i \ -*> d , \ Or . d 2 Wb , 

I j,k,s 



S 



f dv j d2w k „ r f dm dvj dw k 

= L ?>W k e^ v = -J D 

This proves the invariance property ([X, Y - ] |Z) = [X\ [Y, Z]) in case (iii). 

Now we are going to prove the following 
Theorem 5. The generalized Euler equation 

X t = -[AX,X] (4.3) 

for the Lie algebra g(r) with the inertia operator A given by (3.9) yields the equations of 
magnetic hydrodynamics with asymmetric stress tensor (2.2). 

Proof. We write X = ^] Bi-Jj^- + ^2 vjdxj. We will fix representatives of classes of 1-forms 

i 1 j 

modulo dQ°(D) by imposing a condition ^ = 0. Then we have 

.i 3 

and 

^ l dxidxj ^ 3 dxj dxi ^ dxj dxidxk 

Note that the terms Bid(Bi) = ^d(Bf) and Vid(vi) = \d{vf) are full differentials and thus 
may be dropped. 

Substituting the obtained expression into the generalized Euler equation (4.3) and 
collecting terms at dxj, and taking into account that equality of 1-forms is taken 

modulo dQP(D), we get the following system of PDEs: 

dvj ^ dvj ^ dBj ^ dBi d 2 v k dp 



dt ^ 1 dxi ^ 1 dxi ^dxkdxidxj dx^ 

i i i,k J J 

dB j = Y(^ dBj B % dVj 

dt ^ \ 1 dxi 1 dx l 
12 



Rewriting this system in a vector form with the vector fields we get (2.2). 

Finally let us prove Theorem 3 and establish the energy and the cross-helicity con- 
servation laws for MHD equations with the stress tensor (2.2). We will in fact obtain 
Theorem 3 corollary of the following general 

Theorem 6. Let q be a Lie algebra with a non-degenerate symmetric invariant 
bilinear form (-|-). Let A be an involution of q preserving the invariant form, 

A: g^g, A 2 = Id, (AX\AY) = (X\Y) for all X,Y eg. 

Then the generalized Euler equation X t = — [AX, X] has the following two first integrals: 

(AX\X) = Const (4.4) 

and 

(X\X) = Const. (4.5) 

Proof. Let us evaluate (AX\X): 

^(AX\X) = (AX t \X) + (AX\X t ). 

Taking into account that (X\Y) = (AX\AY) and A 2 = Id, we get that (AX t \X) = 
(A 2 X t \AX) = (X t \AX). Thus 

^(AX\X)=2(AX\X t ). 

Substituting the right hand side of the generalized Euler equation for X t we obtain 

^(AX\X) = -2(AX\[AX,X]). 

By invariance of the form we get 

(AX\ [AX, X]) = ([AX, AX] \X) = 0. 

Thus ^ (AX\X) = and (4.4) is established. 

The second conservation law (4.5) is obtained in a similar way: 

- (X\X) = 2 (X t \X) = -2 ([AX, X] \X) = -2 (AX\ [X,X]) = 0. 

This completes the proof of Theorem 6. 

Note that for (4.5) we may drop the requirements that A preserves the invariant form 
and is an involution. 
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We obtain Theorem 3 as an immediate corollary to the previous Theorem, noting that 
the bilinear form (3.7) on g(r) is invariant by Proposition 4 and the inertia operator (3.9) 
is an involution and preserves this form. 

We can see that for X = ^-Bigfr + ^2 v j^ x j^ the & rs ^ integral (4.4) becomes the 

* 1 3 

energy conservation law: 

(AX\X)= I ^v t {x) 2 + J2 B i(*) 2 dV = Const, 

J D i i 

and (4.5) becomes the cross-helicity conservation: 

{X\X) = 2 [ y2vi(x)Bi(x)dV = Const. 
Jd . 
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